main.R

#
# COMPARA DATOS DE MOVILIDAD CON LAS MUERTES POR COVID EN LA CIUDAD DE BS AS
#
# 2021-04-19  version actual:   2021-05-04 
# GRUPO:  
#       Cristian Veiga
#       Cristian De Blasis
#       Gabriel Guiño
#       Guillermo Millich
#       Nicolas Gentile
################################################################################

source("check_packages.R")
check_packages()

## [1] TRUE

city = "Buenos Aires"
trans = "walking"

#Obtiene dataframe de movilidad
# 'dynamic' lo toma desde la pagina
# 'manual' lo toma desde el archivo en el raiz con nombre "applemobilitytrends.csv"
source("mobility_df.R")

## Warning: package 'tidyr' was built under R version 3.6.3

movilidad<-get_mobility_df(city,trans,'dynamic') 

## https://covid19-static.cdn-apple.com/covid19-mobility-data/2107HotfixDev16/v3/en-us/applemobilitytrends-2021-05-09.csv

## Warning

## VERSION DE LA MUESTRA DESACTUALIZADA  17

## https://covid19-static.cdn-apple.com/covid19-mobility-data/2107HotfixDev17/v3/en-us/applemobilitytrends-2021-05-09.csv

#Obtiene dataframe de COVID death
source("covid_df.R")
COVID_19<-get_covid_death_per_day()

## Warning: package 'readr' was built under R version 3.6.3

## Warning: package 'dplyr' was built under R version 3.6.3

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

#Unir los dos dataframes, se usa la fecha como Key (tienen el mismo nombre)
final<- merge(movilidad, COVID_19)

# Grafica las muertes vs la movilidad
source("grafica.R")

## Warning: package 'ggplot2' was built under R version 3.6.3

## Warning: package 'patchwork' was built under R version 3.6.3

## Warning: package 'plotly' was built under R version 3.6.3

## 
## Attaching package: 'plotly'

## The following object is masked from 'package:ggplot2':
## 
##     last_plot

## The following object is masked from 'package:stats':
## 
##     filter

## The following object is masked from 'package:graphics':
## 
##     layout

graficar(final)

# Se observa una reducción desde el principio de Abril donde coincide con los
# siguientes hitos:
# 
# 1/4 al 4/4 Semana Santa (caja roja)
# 8/4: Decreto de reducción de movilidad hasta las 23hs (linea punteada)
# 15/4: DNU de reduccion de movilidad 20hs y clases virtuales, etc (linea punteada)

# Dado a la inrcia del sistema las muertes siguieron aumentando, por tal razon 
# se realizará otra gráfica con las muertes retrasada 15 días.

#Se genera una funcion que genere un retraso en la serie
lag_serie <-function(serie,lag){ 
    x1<-serie[lag:length(serie)]
    # Dado a que la funcion vector genera el valor 0 como default, se hace
    # su propia división para que aparezcan valores NAN correspondientes a las 
    # muestras futuras
    
    x2<-vector(mode = "numeric", length = lag-1)/vector(mode = "numeric", length = lag-1)
    return(c(x1,x2))
  }

lag=15
# Se realiza la grafica con un delay de 15 días
final$MUERTES_POR_DIA=lag_serie(final$MUERTES_POR_DIA,15)
graficar(final)

# Para confirmar un quiebre de tendencia se deberá recabar información de los 
# días subsiguientes a fin de poder minimizar erroes de diagnostico. 
#A priori, se visualiza una disminucion signifiativa en la muertes
# por día.